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Abstract 



We calculate finite mass corrections to leptonic decay constants in HQET 
using QCD sum rules. The results are given to two-loop accuracy includ- 
ing renormalization group analysis. The contribution of unknown two-loop 
anomalous dimension matrix-elements is estimated. We find that the l/mg 
corrections to fp, the decay constant of a pseudoscalar meson, amount 
— (23 lb 4 lb 2)% for the B meson. The first error comes from the sum rule 
uncertainty, the second one from unknown or not accurately known matrix 
elements in HQET. We obtain a ratio of vector to pseudoscalar decay constant 
of 1.00 lb 0.07 + 0{l/ml)) for B mesons. 
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I. INTRODUCTION 



QCD in the limit of infinite heavy quarks has received much interest over the past years 
(cf. and references therein). The theory simphfies considerably in that limit and though 
it still does not allow a complete calculation of hadronic matrix elements, it poses some 
conditions on and gives relations between them that can be used to reduce the model- 
dependence of the analysis of experimental data. Its phenomenological relevance is however 
limited by the fact that quark masses are finite. Generally, matrix elements determined in 
the heavy quark effective theory (HQET) receive corrections due to the finite quark mass. 
Since we are dealing with a field theory that has to be renormalized, the corrections are not 
pure powers, but are modified by large logarithms of type In rriQ / fi where mq is the heavy 
quark mass and fi some hadronic scale. Physically, the logarithms are due to the radiation 
of high energetic gluons off the heavy quark. Although in certain cases some terms in the 
expansion can be shown to be small or absent (as is the case with Luke's theorem 
in general they will introduce new unknown quantities order by order in the heavy quark 
expansion. If HQET is supposed to yield results of phenomenological relevance, it is thus of 
paramount importance to keep those corrections under control. 

In the present paper we address the question of l/mq corrections to one specific matrix 
element, the leptonic decay constant /m of a pseudoscalar meson M with momentum p, 
built of the heavy quark Q and the light anti-quark q: 

(0|g7M75Q|M(p))=^/MP^. (1.1) 

From the analysis of potential models [0], it was known long since that F = fMy/fnq 
approaches a constant as itlq becomes infinitely heavy. The next step in unravelling the 
mysteries of J'm was the discovery that the simple square-root like behaviour gets modified 
by logarithmic terms @]. After the development of the formal framework of heavy quark 
QCD as an effective theory 0, F was soon subject to numerical studies on the lattice 
and by QCD sum rules [0-0. After some fluctuations, the results now agree within the 
errors: 

F(/i = ms) ^ (0.4 - 0.6) GeV^/l (1.2) 

The value of /s, the decay constant of the B meson, that can be derived from the above 
value is however considerably bigger than values obtained in "full QCD" without reference 
to HQET [H,|lOl- Since fs is an important quantity in many applications like Bq-Bq mixing, 
e.g. (cf. [0), its knowlegde is of great importance and thus requires a careful analysis of 
finite mass corrections to F. According to estimates of various authors, l/mq corrections 
to the heavy quark limit amount to about 20% 0,00 • 



In a previous paper |0] we have treated higher order corrections based on a "higher twist" 
expansion of a QCD sum rule for Jb- An explicit expansion in inverse powers of ttlq did 
not prove feasible since the expansion of all parameters inherent in the sum rule analysis in 
inverse powers of rriQ is by no means unique. Instead we defined the "correction term" to 
be the difference between QCD and HQET result. In the range mq > 3 GeV the corrections 
proved to be dominated by the linear term in the heavy quark expansion. On the other hand, 
in an attempt was made to compute 1/mq directly in HQET, but still relied heavily on 
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the expansion of sum rule parameters in l/mq. The resulting l/mq terms are about a factor 
of three larger than the correction obtained in . Thus we feel the necessity to complement 
our old calculation by a corresponding one where the l/mq terms are investigated entirely 
in HQET. Moreover, such an calculation offers the possibility to include next-to-leading 
order renormalization group effects studied recently in |T^. 

Our paper is organized as follows: in Sec. II we recall the formal framework of HQET 
and give necessary definitions, Sec. Ill deals with the heavy quark expansion of correlation 
functions, in Sec. IV we attempt an approximate non-perturbative calculation of the rele- 
vant matrix elements in HQET by QCD sum rules. The last section, Sec. V, is devoted to 
the presentation of the results and their discussion. 



II. PRELIMINARIES AND DEFINITIONS 

Before starting with "heavy" calculations, let us recall the main features of HQET. The 
heavy quark field hy is defined by 

P+Q{x) = e-™«(^-^)/i„(x), (2.1) 

where 

P+ = l(l + |^) (2.2) 

is the projector on the upper components of the heavy quark field Q{x), mq is the mass of 
the heavy quark, and v its four-velocity. The contribution of the lower components to the 
QCD Lagrangian can be integrated out and the new effective Lagrangian be expanded in 
inverse powers of the heavy quark mas^ |n,15]: 



Chqet = Ki{v-D)K + /C + C^^^{^i)S + 0{l/ml). (2.3) 

Zmq Zmq 

At order l/mq, there are two operators]^ the kinetic energy operator 

/C = K{iDfK (2.4) 

whose matrix element over hadron states with equal velocities is related to the kinetic energy 
of the heavy quark within these states. Due to reparametrization invariance K, does not get 



renormalized [16|, whereas the chromomagnetic interaction operator. 



S = h.<y^,gF^'K (2.5) 



^By the heavy quark mass mq we understand the pole mass whose precise definition will be given 
in (iH). 



third one, hy[iv-D)'^hy, can be shown to be of higher order in l/mq by virtue of the equation 
of motion of the heavy quark field. 
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acquires an anomalous dimension 'y™'^^ = 60^/ (4vr) +0{a1) that gives rise to the Wilson 
coefficient Cniag(/i)- 

In HQET, the heavy quark does not take part in the dynamics of a process that is 
essentially determined by the hadron's light degrees of freedom. As a consequence, the 
heavy quark spin becomes a "good" quantum number and the theory exhibits the famous 
spin-symmetry0 which leads to a factorization of matrix elements into terms describing 
the spin-structure of the involved hadrons (and their interaction) and terms describing the 
dynamics of the process that is independent of the heavy quarks. The formal tool to handle 
that factorization is the so-called trace-formalism |]T^. In this formalism and using the 
notation of 0, we define, making explicit the dependence on the renormalization scale /j, 
whenever necessary, 

{0\qTK\M) = ^F{fi)TT{TM}. (2.6) 
M. is the spin wave function of the HQET meson state | M ) : 

A ^ / r, i ""^75 for = 0~, ,^ 

M = V^P^[ ; forJ^ = l-: ^2.7) 

Here is the polarization vector of the vector meson. 

At order l/mq and neglecting radiative corrections, the leptonic decay constant Jm is 
determined by (cf. 0) 

/m =f(i- cLmJ^ + ^ + dM—\ (2.8) 

A is defined as the mass-difference of meson mass tum and quark mass rriQ in the heavy 
quark limit, niM = fnq + A + ©(l/mg). Gk and G-e are defined by 

{Q\ijd^xlC{x) g(0)r/i^(0) |M) = ^ F(/i)Gx(Ai)Tr{rM}, 

(0|i|ci^a;5(x)g(0)r/i„(0) |M) = ^ F(/x)Gs(/^)2rfMTr{rA^}. (2.9) 

As already mentioned, at order l/mg the spin-symmetry breaks down and matrix elements 
become dependent on the spin of the heavy particle as characterized by a scalar dM defined 
by 

P+a^fsMcr''^ = 2dMM, (2.10) 

yielding 



^Of course the symmetry breaks down at order l/mg, but in a so to speak controlled manner 
that is induced by S and only affects the traces over spin wave functions. 
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- { -1 for r = 1-. ^^-^^^ 

Eq. (|2.8| ) contains all information about Jm to order l/rriQ and to lowest order in the 
strong coupling, but involves the three new matrix elements A, Gk and Gj^. For A there 
exists a rigorous lower bound of 237 MeV [0 that is however questioned in |T9[. Neglecting 
l/mg effects, the values of the b-quark mass quoted in literature, = (4.6 — 4.8) GeV 
yield A = (0.5 — 0.7) GeV. More refined analyses in HQET tend to the somewhat 



smaller value A = (0.5 — 0.6) GeV which we will use in our paper. 

The available values of Gk and Gs rely exclusively on the analysis given in |^ that is 
similar in spirit to the one we are going to undertake but quite different in detail. We will 
come back to that point later. The results of read (at a renormalization scale of 5 GeV): 

Gk ^ -4.4 GeV, Gs ^ -0.05 GeV. (2.12) 



III. THE HEAVY QUARK EXPANSION OF CORRELATION FUNCTIONS 

Our calculation of Gk and G^ relies on the analysis of two-point functions of the heavy- 
light current Jr = qTh^ with insertions of /C and S, respectively. Since we will include radia- 
tive corrections to 0{as), an investigation of the renormalizat ion-group behaviour proves 
necessary. In addition, at the required level of accuracy we need to match the effective 
theory to QCD, i.e. we have to require the equality of matrix elements in HQET to the cor- 
responding ones in QCD at the matching scale /i = mq. As a by-product of this procedure, 
we will be able to check the heavy quark expansion explicitly at the two-loop level at the 
example of the two-point function of pseudoscalar currents. This will be done in the next 
section. 

Let us first introduce a generic notation for correlation functions of two composite oper- 
ators Oi and O2, either containing a heavy quark field, 

{O1O2) =ijd^xe''^'''' {0\TOi{x)Ol{0)\0), (3.1) 

with four-momentum ujv, v being the velocity of the heavy quark. Analogously, we define 
for the zero-momentum insertion of the operator X: 

{O1XO2) =i^Jd^xd^ye"'''<^-y^0\TOi{x)X{0)Ol{y)\0). (3.2) 

Employing the trace-formalism for heavy-light currents Jr^ = qTihy we can writeQ 

{JrM = ln{u;,ft)TT{T,P+T2^}. (3.3) 



^Our definition deviates formally from the conventional one ( Jr^Jpa) = — ^ n(a;) Tr(riP-|_r2), 
but coincides for "spin eigenstates" like Fj = ^75 or Fj = 7^ — with Fj|^ = — |^Fj and can be 
directly inferred from the structure of the relevant loop-integrals. 
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Recall that 11(0;, /i) is independent of both Fi and but depends on the renormalization 
scale /i due to the non-vanishing anomalous dimension of Jp. With X = K, or S, (|3.2| ) can 
be expressed as 

( Jr,/C Jr, ) = Tk{uj, fi) Tr{riP+r2^}, 

(Jr.SJr,) = T^iuj,fi)dMTr{T,P+T2ij}. (3.4) 

In QCD, the renormalization-group invariant two-point function of pseudoscalar currents 
with momentum q is defined as 

Hs = tld^xe'''^0\Tm^q{x)tTM^)my^QiO)tl,qiO)\0). (3.5) 

Here the running MS mass m-^ enters whose relation to the pole mass is given by 



l + ^P-21n^ +0f 

TT \ 3 fi J V 



(3.6) 



lis was calculated to two-loop accurcay in the deep Euclidean region <^ in |]2T]] for 
arbitrary quark masses. Once the heavy quark expansion of the pseudoscalar current is 
known, the expression for lis can serve as check for both the cancellation of fi dependent 
terms in the effective theory and the values of the matching coefficients, i.e. the values of 
the short-distance expansion coefficients at yU = rriQ. 

As for the pseudoscalar current, its expansion reads (following the notation of ||13||): 



/ 1 2 1 2 \ 
^MS^^75Q = niQ C(^)J,^3 + ^ + ^ Y.Ml^)Tk] ■ (3.7) 

\ Q i Q k / 

The above formula needs some explanation. First, "=" means equality only after taking 
matrix elements or insertion in Greens functions. Next, we have kept the multiplicative 
mass term in order to indicate that it affects the coefficient functions on the right hand side 
by virtue of Eq. (|3l6|). The leading order operator in HQET is the current Ji^^; at order 
l/mg, there are two local operators, 

Oi = qi'j5iphv, O2 = q{-iv-'D)i'y5h^, (3.8) 
and two non-local ones due to the insertion of the Lagrangian into the leading term: 

Ti = iJd^xTJi^^{x)}C{0), T2 = iJd'^xTJi^^{x)S{0). (3.9) 

The operators have non-vanishing anomalous dimensions and mix under renormaliza- 
tion. The mixing is described by an 4 x 4 matrix 7^^ and the coefficient functions 
{Bi, B2, Ai, A2} = C obey the renormalization-group equation 

(,A_,Ps\^.0. (3,10) 
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Since /C does not get renormalized, it follows Ai{fi) = C{fi). For A2, we have A2{fi) 



C(A*)Cmag(Ai) due to the non-zero anomalous dimension of S (cf. (|2.3| )). Cmag "was calculated 
Mm to be 



m 



= ! + -(- + - In — 



(3.11) 



In order to determine the remaining coefficients, we proceed like Ref. [|I^ and calculate 
matrix-elements between on-shell quark states. The corresponding diagram is shown in 
Fig. |1|. The on-shell spinors obey the Dirac equations -^uq^p) = mqUqij)) and fuq{p') = 0. 
The heavy quark momentum p can be written asp = rnqv + k where the off-shell momentum 
k obeys 2mQV-k + k^ = in order to keep uq on-shell. For the on-shell spinor in HQET the 
Dirac equation reads '^Uh = 0, and it follows the relation uq = {1 + ^/ {2mq)}uh. Taking into 
account only terms constant or linear in p', loop-integrals in HQET vanish in dimensional 
regularization since there is no massive scale they could depend on. The diagram in Fig. |l|, 
wave function renormalization and terms coming from the replacement of the running MS 
mass by the pole mass yieldQ 




(3.12) 



The expression for C(/i) is well known and can be found in a somewhat different form in ||23 
apart from a slight error that was corrected in [|^,^ . The first equality in ( p.l2| ) is consistent 
with the requirement of reparametrization invariance, the expression for B2 is new. 
The expansion of ( p.5|) to order 1/mq is now straightforward to derive. Using 



(3.13) 



( -h^iDjvi ) = ^^n(^) Tr{(7^ - t;„|^)riP+r2}, 
o 



the insertion of (13.71) into 



yields 



{C'ifi) - B2{fi)C{fi)}ujU{uj, /i) + C'{fi)TK{uj, fi) 



+3C2(/i)Cmag(/i)7s(t^, /i) 



(3.14) 



This expression is valid up to constants in u that are related to contact terms and cannot 
be expanded in 1/mq. The //-dependent terms on the right hand side must cancel in order 
to keep the result renormalization-group invariant. This cancellation will be checked in the 
next section. 

One could now proceed by solving the renormalization-group equation (|3.10|) as was done 
in 1^. For our purposes, however, it proves more appropriate to investigate the behaviour of 
the two-point functions V = {uII{u),Tk{u),Ty;{uj)} under renormalization-group. V obeys 
the renormalization group equation 



^We use anticommuting 75. 
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f^^ + l] V(/i) = 0. 



(3.15) 



The 3x3 anomalous dimension matrix 7 can be obtained from the renormahzation constant 
Z that relates the bare Vbare to the renormahzed one, Vbare = ZV. Restricting ourselves to 
leading logarithmic accuracy (superscript LL), 7^^ can be obtained as 7^^ = lime^o2eZ in 
the MS scheme^] using dimensional regularization in D = 4 + 2e dimensions. For ^ 0, V 
can be treated within perturbation theory and we find 



/ -8 



7 



LL 



as 
Air 



V 



32 Q 
9 ^ 





-2; 



(3.16) 



To that accuracy, the solution of Eq. ( ^.15| ) is 

V(mQ) = exp 



Wo^""") ^^^^ 



(3.17) 



with X = as{iJi) / Oisirnq). [3 = —Qs [/3otts/(47r) + (3i{as/ {A'n)Y + 0{a1)] is the [3 function of 
QCD. The matrix 7 cannot be diagonalized, but transformed to its Jordan form 7j via 



7j = T ^^T 



with 




7j 




Using 



exp 



a 1 
a 



exp(aj 



1 1 
1 



we find 



exp 



7 



In a; 



"3/3o 




1 



\ 




V -if (x-3//3o -1)0 



X 



-3/A) 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



As expected, Vi renormalizes like n(ci;). For V2 and V3 we can infer the renormalizat ion- 
group invariant quantities 



But note that the result is renormalization-scheme independent to that accuracy. 
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In order to extend the analysis to next-to-leading order (superscript NLO), one needs to 
include the anomalous dimension matrix 7 to two-loop accurcay. For Vi one finds []7|,P|,pi 

= fl - ^HL ^] Vi(/i), (3.23) 



27r 



where 

^HL / HL 



^HL = ^ ?k-^ • (3.24) 



2/3o V75'^ A 



0, 



^HL _ ^HL^^^j^^^-j _|_ ^HLj^^^j^^^^ j2 _|_(^(^q,3-j jg auomalous dimension of the heavy-light 
current Jr calculated in p3| , ^ . As for V2 and V3, some of the involved elements of 7 are 
still unknown to date. According to |13| one obtains 

% = [a.(/i)]^/^o fl - 5„L ^] Uif^) - V^i(/i) ( ^ ln«.(/i) + |t34 "^^^^^ ^ 



3/5o ■ 3 - 



V, = i„,(,)i"«> (1 - 2s„. _ s„.. ^^) ^3(,<) - ^ v.(.) (1 ^ 



(3.25) 



S'mag is the analogue of S'hl were 7 has to be replaced by the anomalous dimension 7^^^^ 
of the operator S. Unknown are the matrix elements T32 and T34 as well as 7™^^. 

We now want to relate the two-point functions we have been dealing with in this section 
to the l/mq corrections to F we are interested in. To that end, we saturate the correlation 
function lis by hadronic states. Making explicit only the insertion of the ground-state, the 
meson M, it follows with ([TTI) 

^ m|(m|^-g2) 

where the dots denote insertions of radial excitations and multi-particle states. Analogously, 
we find from ( p.9|) , choosing /i = mg: 

r,(.,™«) = £!i^^2i^ + ... (3.27) 

[A - to) 
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For Tk and the dots contain also additional double poles at = A that are related to 



l/mQ corrections to A and are treated in detail in |26 



Identifying the leading terms in l/mg, we recover from ( p.l4| ) the well known relation 

(3.28) 



for the pseudoscalar decay constant fp p3| 
At order l/mq, we find 



fpy/m^ = F{mQ) 1 - 



2 a 



■{jnQ)\ 



1 + 



TT 



+ 



3/ ^13a>^\ 



6 



TT 



(3.29) 



From an analogous consideration of the vector current, it follows for the vector decay con- 
stant fv imi 



TT 



TT 



A + -G'j^(mQ) 



- 1 



13 as("^Q) \ 

6 TT i 



G's("^Q) 



(3.30) 



These expressions are the generalization of (|2.8| ) to order a^. Eq. (|3.29|) coincides with the 
corresponding expression derived in |T^. 

The ratio of vector to pseudoscalar decay constant is then given by 

!v^^_2ajjr^^^l^2 1^'. (3.31) 

fp 37r niQ \ -n ) \ 2 -n ) ] 



From ( |3.25| ) we now can construct the invariant quantities 
Gk = GM-K 



'4 lna.(^) + ^T34^' 

3/5o 3 TT 



[«,(/i)]-3/'3« 1-5, 



mag 



47r 



27 



TT 



(3.32) 



IV. SUM RULES FOR Gk AND Gs 

We now turn to the main part of this paper: the derivation of QCD sum rules for Gk and 
Gy,. As mentioned in Sec. I, a calculation of Gk and from first principles is not feasible 
to date, nor is it for hadronic matrix elements in general. Instead one has to rely either on 
models invented ad hoc whose relation to the underlying field theory remains unclear or one is 
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restricted to numerical approximations to QCD that allow the inclusion of non-perturbative 
effects in a systematic way. One of these appoximations is lattice QCD whose results can in 
principle be refined to arbitrary accuracy, but in practice are limited by available computer 
power and time. An other, less ambitious method is provided by QCD sum rules (cf. ||27 



for reviews). Originally invented for the calculation of vacuum-to-meson amplitudes in 
QCD pSf, it found application likewise to meson-to-meson transition amplitudes ||29| and 



could also be adapted to the study of HQET |]7|,P|pi|. QCD sum rules rely on the field- 
theoretical aspects and features of QCD and were designed to make maximum use of known 
manifestations of non-perturbative QCD. Their advantage over lattice calculations is mainly 
their simplicity, whereas the main disadvantage may be seen in the fact that a systematic 
refinement to a given accuracy is not possible. The best accuracy one can hope to achieve 
is at the (10 - 20)% level. 

Without going too much into details, we want to present here only the gross features 
of the method. QCD sum rules are based on the operator product expansion (OPE) of 
some correlation function, say n(co'), in terms of the so-called condensates, vacuum expec- 
tation values of certain operators that vanish when taken over the perturbative vacuum, 
but acquire non-zero values in the physical QCD vacuum and thus characterize its non- 
perturbative nature. The expansion is done in the not so deep Euclidean region of uj where 
non-perturbative effects are already noticeable, but do not dominate yet: 

En('^)(^)(a) =^n(")(^)(o|a|o) = n(a;). (4.i) 

n n 

Examples for the vacuum condensates (On) entering the OPE are (Oi) = U, the unity 
operator, {O3) = {qq), the (light) quark condensate, {O4) = {asC^/n), the gluon condensate, 
and (O5) = {qgsCr^^G^'^q) , the mixed condensate. The coefficient functions n'^"^(a;), also 
called Wilson coefficients, can be calculated within perturbation theory and are determined 
by the short-distance behaviour of Il{uj). 

On the other hand, 11 can be expressed via Cauchy's theorem in terms of its singularities, 
which are situated on the positive real axis in u and are due to bound states causing poles and 
multi-particle states causing cuts. If we use as an input the information, that there is a pole 
due to the lowest-lying one-particle state, and its position, we can fit its residue, the square 
of the vacuum-ground-state transition amplitude, from the knowledge of 11 in the not so 
deep Euclidean region. The short-comings, however, of such an approach are obvious: first 
we have to take it for granted that the OPE exists and makes sense. Besides, our knowledge 
of n in the not so deep Euclidean region is not complete since one has to truncate the OPE 
at some point, which leads to errors in the determination of the residue. Moreover, there is 
not only the one-particle ground state, but also radial excitations and multi-particle states 
contribute to the singularities of 11. Although it seems natural to assume that contributions 
of multi-particle states be suppressed by creation amplitudes, one still has to deal with 
both the unwanted contributions from radial excitations and unknown terms in the OPE 
stemming from higher-dimensional condensates. 

To this end, it proved extremely useful in many applications to subject 11 to a Borel 
transformation, a tool, that likewise was introduced in [|28|. The Borel transform BF of a 
function F depending on the Euclidean momentum ue = —oj is defined as 
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1 

BF= lim ^-^^ F. (4.2) 

u}E/N=t fixed 

The dimensionful quantity replacing uj is t, the Borel parameter. Apphed to a dispersion 
relation over the physical spectral function pP^^^^, e.g., the Borel transformation yields 

B jds = \ Jdsp^'^^'is) e-'l\ (4.3) 

Thus contributions of p^'^^^^s) become exponentially suppressed for large s, the suppression 
being controlled by the value of t. Contributions of higher dimensional condensates to 
the OPE (|4.1|) get suppression factors in powers of 1/t and factorials. Roughly speaking, 



the Borel parameter serves as weight between perturbative (oc (0i)) and non-perturbative 
contributions. 

Yet, one has to model pP'^^*^(s) for large values of s. At that point the quark hadron 
duality hypothesis enters. It states that integrals over pP'^y'^ can be approximated to good 
accuracy by integrals over the perturbative spectral function pP'^''*, multiplying (Oi), to wit 

Jds pP^^^s) e-'/' ^ JdspP'''\s)e-'/'. (4.4) 

Thus a suitable model for pP^y*^ for large values of s is pP^^'^(s) = pP®''*(s)0(s — c^o) where 
integration is understood afterwards, uq is the so-called continuum threshold and is un- 
fortunately not related to fundamental quantities of QCD (or HQET), but should be less 
or equal to the position of the lowest-lying excitation, i.e. c^o ~ 1 GeV. In the numerical 
analysis of sum rules one thus has to look for a value of t^o and a sufficiently large range of t 
where the sum rules becomes stable in these parameters. The accuracy of about (10 — 20)% 
mentioned above can, however, be worsened by a possibly strong dependence of the result on 
the precise values of the Borel parameter and the continuum threshold. Stated differently, 
too strong a dependence on that parameters is an indication for the unreliabiliy of a sum 
rule. As we will see in the next section, that is fortunately not the case with the quantities 
we are interested in. 

Keeping all these caveats in mind, we proceed with the derivation of QCD sum rules for 
Gk and Gj^. We write 

T.^fd.eiM, T,= [ds>^, (4.5) 

J S — UJ J S — 00 



do the OPE at the level of spectral densities and thus define 

PK = pT + P^Sm + pS^ /- + pf{qgsa,.G^''q) + • • • (4.6) 



and correspondingly for pg. We use Feynman gauge for all diagrams involving only pertur- 
bative gluons; for the calculation of non-perturbative diagrams we choose the coordinate 
gauge [|30|. The sign of the strong coupling is defined by = — We ap- 



ply dimensional regularization in D = 4 + 2e dimensions and use the MS renormalization 
scheme. For the perturbative contribution, we take into account the diagrams shown in Fig. 
H and find (at the renormalization scale p): 
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TT \ 6 3 P J ) 



9 3 ^ 



TT 



TT 



P 



e(s). 



(4.7) 



The loop-integrals involved in the calculation can be found in p6[, e.g. The diagrams for 
the contribution of the quark condensate are shown in Fig. ^ and yield 



(3) 

Pk 



0, P§' 



2 



e(s). 



IT 



(4.8) 



In addition we have taken into account the contributions of the gluon and the mixed con- 
densate shown in Fig. ^ that read 



(4) 

Pk 



(4) 



-48 



(5) 

Pk 



32 



48 ' ' 48 
From the result for 11 (cj) given in |l7|J^,p^ we obtain 



(4.9) 



PojU 



1 + 



asifi) /17 4 



TT 



— + -vr" 
3 9 



21n^j|e(s) 



Y^S{s){qg,a^,G^''q). (4.10) 



We can now check all the formulae obtained so far by means of Eq. ( |3.14 ). The expression 
for lis including non-perturbative corrections can be found in [T^ , e.g. Here we concentrate 
on the perturbative part calculated in |21|. For its expansion to 0(l/mQ), we recall that 
( p.l4| ) is valid up to constants, so we make frequently use of partial integration to obtain 
the spectral density. Expanding the four-momentum as = tUq + 2mQUj + cj^, we findQ 



pert 



Q 



a 
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2s 



27r2 [ 
1 



+ 



1 + — 2 In + -7r 

TT \ 3 irtQ 9 

2uj 



rriQ TT" 



3 as 
+ V 



-14 + 3 In- 



TT 



e{s). 



(4.11) 



Inserting C(/i), -Bi(/i) and B2{n) from ( p.l2| ) and pj^n*; pW^ P?'^^ from above into the 
right hand side of (^.141), we indeed recover (|4.11|). 



^For the expansion of the full lis, note that the radiative correction to the Wilson coefficient 
multiplying the quark condensate was first calculated in and reads before application of the 
Borel transformation: 

nf^(^g) = (mj^(gg)) 
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There remains only one slight complication to cope with, that is the structure of Tk 
and Ts expressed in terms of hadronic matrix elements. As mentioned at the end of the 
last section, both quantities develop apart from the single poles at = A we are interested 
in also double poles that stem from the contribution of three-point functions^ and can be 
obtained by inserting hadronic states between each of the operators in (Jr^OJrj)- With 
K = (M|/C|M)/(2mj\//) the invariant quantity Tk defined in (p.25|) can be expressed as 



F^Gkt F'^K 

We can however isolate the single pole contribution by first applying a Borel transformation, 

BTk = ^ F'ke-^/' + i F'Gkc-^/' + ... 

= ^kBfl + GKBfl+ (4.13) 

where in the step from the first to the second equation we have made use of the Borel 
transform of Eq. (|3.27|) . Taking then the derivative with respect to t yields 



^ d IBTk 

GK = ——r^. 4.14 
dt EH ^ ' 



The corresponding equation for Gs is 



^ 2dt BH ^ ' 



The residues of the double pole contributions were calculated in ||2^ . Eqs. (|4.14|) and ([4.1 5|) 
are the final sum rules we are going to evaluate in the next section. 

Let us close with a few remarks about the calculation performed in 0. In our approach, 
the contributions of 0(1) are completely separated from the ones of 0(l/mQ), because the 
matrix elements considered are different. So from our point of view it is not allowed to 
introduce a shift of the Borel parameter of order l/mq in the sum rule for Jm as it was 
done in ^ since this would mix different orders in the inverse heavy quark mass. Since 
we have derived sum rules for Gk and Gs separately, there is no need to shift the value of 
the continuum threshold, either. Furthermore, the double pole contributions escaped the 
attention of Ref. P and thus the results obtained there are not values for Gk and G^, but 
mixtures of Gk and K and of Gs and {M\S\M). 

V. RESULTS AND DISCUSSION 



In the numerical evaluation of the sum rules ( |4.14| ) and ( [1.15| ) we use the following values 



of the condensates at the renormalization point = 1 GeV 



®Such a mixing of two- and three-point functions was first considered in [31 
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{qq) = (-0.24 GeV)^ {^G^'j = 0.012 GeV^ 

{^gsO^.G^'q) = 0.8 GeV^ {qq). (5.1) 
For scales different from 1 GeV, we scale the condensates accordingly. For a^, we use 



Air 



2/3oln(/i/Ag|)) 



^ ln(21n(/./A^0) 
/5o^ 21n(/./Ag|)) 



(5.2) 



with = 260 MeV, corresponding to a,(5 GeV) = 0.20. 



MS 



In Fig. 1^ we show Gk as function of t for c^o = 1 GeV. Although TV is renormalization 
group invariant to next-to-leading order, we still have to choose one particular scale in 
building this quantity (cf. (|3.25|) ). As argued in [0, the most natural choice is /i = 2t. In 
Fig. 1^ we also show curves for /i = 1 GeV and /i = 2 GeV. The scale-dependence still visible 
indicates the size of higher order effects in a^. also depends on the unknown anomalous 
dimension matrix element T34 whose value is put zero in the figure. Varying the continuum 
threshold within ujq = (1.0 — 1.2) GeV and t within 0.35 GeV < t < 0.6 GeV, as suggested 
by the analysis for F done in |^, and T34 in the range from —1 to 1, we obtain 

= -(1.1 ±0.2 ±0.2) GeV. (5.3) 

The first error indicates the sum rule uncertainty in the "sum rule window" 0.35 GeV < t < 
0.6 GeV, the second one the uncertainty stemming from the variation of T34 which is rather 



large. From ( |3.32| ) we then obtain 

GximQ) = -(1.6 ± 0.2 ± 0.3 ± 0.1) GeV, (5.4) 

where the first two errors have the same meaning as before, the third one stems from varying 
A in the range (0.5 — 0.6) GeV. 

In Fig. 1^, Ge is displayed as function of t for uq = 1 GeV. As for Gk, we have varied /i 
from 1 to 2 GeV. Although here two unknown matrix elements enter, Smag and T32 whose 
values we have put to zero in the figure, their influence on the result is much smaller than 
that of T34 on Gk and we find 

= -(0.19 ±0.06 ±0.02) GeV. (5.5) 

which corresponds to 

G^iruQ) = (0.042 ± 0.034 ± 0.023 ± 0.030) GeV. (5.6) 

The large errors in this result come from a cancellation of the central values of the two terms 
in (PI). 

In the last figure. Fig. ||, we show the resulting l/rriQ corrections for fp^frnp and fy^fmy, 
denoted as dFjF, that are obtained from Eqs. (|3.29| ) and ( p.30| ). They yieldQ 



^Note that the errors of Gj]{mQ) and Gk{t^q) are correlated. 
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fp^ff^ = Fimn) I 1 - - "''^"^^M ( 1 - — (1.07 ± 0.19 ± 0.07 ± 0.02) GeV 
V 3 TT ; V mQ J 




fyJ^ = Fimn) ( 1 - - ^ii!^ ) ( 1 - — (0.78 ± 0.08 ± 0.09 ± 0.04) GeV ) (5.7) 

V 3 TT ; V mg ; 

with the same sequence of errors as before. Note that the error stemming from the un- 
known matrix elements is reduced as compared with ( p.4|) and ( ^.61) and that the sum rule 
uncertainty for /y is smaller than that for fp due to the better stability visible in Fig. ^. 
Finally, we obtain from ( p.31| ) (at the scale rrib) 

^ = 1.00 ± 0.03 ± 0.02 ± 0.02 + 0{l/ml), (5.8) 
Jp 

that is, a nearly cancellation of the ^jvriQ term, indicating that l/'n^g correction terms 
become important. 

The result ( |5.7|) agrees well with our old one obtained from the QCD sum rule for fp [0: 

The remaining slight discrepancy is within the uncertainty of the method and/or could be 
due to 1/?tiq correction terms. Our result also agrees well with values obtained from lattice 
calculations 0. As far as p is concerned, our value for Gxifng) is by nearly a factor three 
smaller, for Gs(?71q) we obtain a different sign. The difference in sign affects mainly the 



ratio fv/ fp, where we find some disagreement with the lattice calculation where fv/ fp 
is found to be 1.11 ±0.05 at the scale of the B meson. This value, however, relies on a linear 
fit through four points in the mass range (1.3 — 2.6) GeV that has to be extrapolated to a 
mass of 5.28 GeV. Fitting a parabel reduces the ratio already to 1.06 instead of 1.11. The 
remaining discrepancy could be due to l/rriQ terms that mimic a l/rriQ behaviour in the 



restricted mass range considered in 



Concluding, we have found the sum rules for l/mq corrections in HQET in complete 
agreement with conventional sum rules. There is however, a discrepancy in fy / fp compared 
to lattice calculations which might be due to ^/ttlq correction terms. 



ACKNOWLEDGMENTS 

It is a pleasure to thank V.M. Braun for useful discussions and collaboration in the early 
stages of this work. This work was supported in part by the German Bundesministerium 
fiir Forschung und Technologie under contract 06 TM 732 and by the GEO science project 
SC1-CT91-0729. 



16 



REFERENCES 



[1] M.B. Wise, in Proceedings of th 6th Lake Louise Winter Institute, Lake Louise, 1991, 
edited by B. Campbell et al. (World Scientific, Singapore, 1991), p. 222; H. Georgi, 
in Proceedings of the TASI 1991, Boulder, 1991, edited by R.K. Ellis, C.T. Hill, and 
J.O. Sykken (World Scientific, Singapore, 1992), p. 589; M. Neubert, Preprint SLAC- 
PUB-6258 (1993). 

[2] M.E. Luke, Phys. Lett. B 252, 447 (1990). 

[3] Ya. I. Azimov, L.L. Frankfurt, and V.A. Khoze, JETP Lett. 24, 338 (1976); E.V. 

Shuryak, Nucl. Phys. B198, 83 (1982). 
[4] M.A. Shifman and M.B. Voloshin, Sov. J. Nucl. Phys. 45, 292 (1987); 47, 511 (1988). 
[5] H. Georgi, Phys. Lett. B 240, 447 (1990). 

[6] A. Duncan et al, in Lattice 92, Nucl. Phys. B (Proc. Suppl.) 30, 433 (1993); C. Alexan- 
drou et al. Preprint PSI-PR-92-27 (1992) (|hep-lat/9211042|) ; C.W. Bernard, J.N. 
Labrenz, and A. Soni, Washington Preprint UW-PT-93-06 (1993) ( |hep-lat/ 930501151 ); 
APE Collaboration (C.R. Allton et al). Preprint LPTENS-93-12 (1993); UKQCD Col- 
laboration (R.M. Baxter et al), Edinburgh Preprint 93/526 (1993) ( |hep-lat/93(MM]| ). 

[7] E. Bagan et al, Phys. Lett. B 278, 457 (1992). 

[8] M. Neubert, Phys. Rev. D 45, 2451 (1992). 

[9] M. Neubert, Phys. Rev. D 46, 1076 (1992). 
[10] For a recent overview over QCD sum rule calculations, see C.A. Dominguez, Talk given 
at the Third Workshop on the Tau-Charm Factory, Marbella, Spain, 1-6 June 1993, 
Preprint SISSA Ref. 139/93/EP (1993) ( |hep-ph/9309260D . 
[11] B-Physics Working Group Report, edited by R. Aleksan and A. Ali, DESY Report 
93-053 (1993). 

[12] V. Eletsky and E. Shuryak, Phys. Lett. B 276, 191 (1992). 

[13] M. Neubert, Preprint SLAC-PUB-6344 (1993) (|hep-ph/9308"369| ). 

[14] A.F. Falk, B. Grinstein, and M.E. Luke, Nucl. Phys. B357, 185 (1991). 

[15] T. Manuel, W. Roberts, and Z. Ryzak, Nucl. Phys. B368, 204 (1992). 

[16] M. Luke and A.V. Manohar, Phys. Lett. B 286, 348 (1992). 

[17] J.D. Bjorken, in Proceedings of Les Rencontres de Physique de la Vallee d'Aoste, La 
Thuile, 1990, edited by M. Greco (Editions Frontieres, Gif-sur-Yvette, 1990), p. 583; A. 
Falk et al, Nucl. Phys. B343, 1 (1990). 

[18] Z. Guralnik and A.V. Manohar, Phys. Lett. B 302, 103 (1993). 

[19] I.I. Bigi and N.G. Urahsev, Preprint CERN-TH.7091/93 (1993) ( |hep-ph/93im7D . 

[20] M.B. Voloshin and Yu.M. Zaitsev, Sov. Phys. Usp. 30, 553 (1987); L.J. Reinders, Phys. 
Rev. D 38, 947 (1988); C.A. Dominguez and N. Paver, Phys. Lett. B 293, 197 (1992); 
P.B. Mackenzie, Talk presented at the Lepton-Photon Symposium, Cornell University, 
Aug. 10-15, 1993, Preprint FERMILAB-CONF-93/343-T (1993) (|hep-ph/9311242| ); S. 
Titard and F.J. Yndurain, Madrid Preprint UM-TH-93-25 (1993) (|hep-ph/931023q ). 

[21] D.J. Broadhurst, Phys. Lett. B 101, 423 (1981). 

[22] E. Eichten and B. Hill, Phys. Lett. B 243, 427 (1990). 

[23] X. Ji and M.J. Musolf, Phys. Lett. B 257, 409 (1991). 

[24] D.J. Broadhurst and A.G. Grozin, Phys. Lett. B 274, 421 (1992). 

[25] D.J. Broadhurst and A.G. Grozin, Phys. Lett. B 267, 105 (1991). 



17 



[26] P. Ball and V.M. Braun, Miinchen Preprint TUM-T31-42/93 (1993) ( |hep-ph/930729TD , 

to appear in Pliys. Rev. D. 
[27] L.J. Reinders, H. Rubinstein, and S. Yazaki, Phys. Rep. 127, 1 (1985); S. Narison, 

QCD Spectral Sum Rules (World Scientific, Singapore, 1989); Vacuum Structure and 

QCD Sum Rules, Current Physics Sources and Comments, Vol. 10, edited by M.A. 

Shifman (North-Holland, Amsterdam, 1992). 
[28] M.A. Shifman, A.I. Vainshtein, and V.I. Zakharov, Nucl. Phys. B147, 385, 448, 519 

(1979). 

[29] B.L. loffe, A.V. Smilga, Phys. Lett. B 114, 353 (1982); Nucl. Phys. B 216, 373 (1983). 
[30] V. Fateev, A. Schwartz, and Y. Tyupkin, Moscow Preprint FIAN-155 (1976); C. Cron- 

str5m, Phys. Lett. B 90, 267 (1980). 
[31] V.M. Braun, A.V. Kolesnichenko, Nucl. Phys. B283, 723 (1987). 
[32] A. Abada et al, Nucl. Phys. B376, 172 (1992). 



18 



FIGURES 




P P' 



FIG. 1. Momentum configuration for tlie matching calculation of the pseudoscalar current, p 
is the momentum of the heavy quark Q, p' that of the massless quark q. Both quarks are on-shell. 




FIG. 2. Perturbative contribution to Tk- Double lines denote heavy quarks, single lines light 
quarks and curly lines gluons. The crosses with circles denote the heavy-light vertex F, the black 
square the insertion of the operator /C. All diagrams not drawn vanish. 




FIG. 3. Diagrams contributing to the Wilson coefficient of the quark condensate in the OPE 
of Tk- Lines with crosses denote vacuum expectation values. 
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FIG. 4. Diagrams contributing to the Wilson coefficient of the gluon and the mixed condensate 
in the OPE of Tk- We use the coordinate gauge with the coordinates specified in the second 
diagram. All diagrams not drawn vanish. 




FIG. 5. Diagrams contributing to T^;. The black square now denotes an insertion of the operator 

S. 
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FIG. 6. Gk as function of t for tuo = 1 GeV and different scales /i. Solid line: /i = 2t, 
long-dashed line: ^ = 1 GeV, short-dashed line: fi = 2 GeV. 
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FIG. 7. Gs as function of t for = 1 GeV and different scales Solid line: fj, = 2t, 
long-dashed line: /u = 1 GeV, short-dashed line: fi = 2 GeV 
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FIG. 8. The l/mq corrections 6F/F to F according to Eqs. ( 3.2g| ) and (|3.30| ) as function of t 



for ujQ = 1 GeV. Solid line: pseudoscalar decay constant, dashed line: vector decay constant. The 
normalization scale is /i = 2t. 
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